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Abstract. Sulfur dioxide (SO 2 ) is a major atmospheric pol- 
lutant with a strong anthropogenic component mostly pro- 
duced by the combustion of fossil fuel and other industrial 
activities. As a precursor of sulfate aerosols that affect cli- 35 
mate, air quality, and human health, this gas needs to be 
monitored on a global scale. Global climate and chemistry 
models including aerosol processes along with their radia- 
tive effects are important tools for climate and air quality re- 
search. Validation of these models against in-situ and satel- 40 
lite measurements is essential to ascertain the credibility of 
these models and to guide model improvements. In this 
study the Goddard Chemistry, Aerosol, Radiation, and Trans- 
port (GOCART) module running on-line inside the Goddard 
Earth Observing System version 5 (GEOS-5) model is used 45 
to simulate aerosol and SO 2 concentrations. Data taken in 
November 2010 over Frostburg, Maryland during an SO 2 
field campaign involving ground instrumentation and aircraft 
are used to evaluate GEOS-5 simulated SO 2 concentrations. 
Preliminary data analysis indicated the model overestimated 50 
surface SO 2 concentration, which motivated the examination 
of mixing processes in the model and the specification of SO 2 
anthropogenic emission rates. As a result of this analysis, a 
revision of anthropogenic emission inventories in GEOS-5 
was implemented, and the vertical placement of SO 2 sources 55 
was updated. Results show that these revisions improve the 
model agreement with observations locally and in regions 
outside the area of this field campaign. In particular, we 
use the ground-based measurements collected by the United 
States Environmental Protection Agency (US EPA) for the 60 
year 2010 to evaluate the revised model simulations over 
North America. 
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1 Introduction 

Sulfur dioxide (SO 2 ) is a trace gas which poses significant 
health threats near the surface, with consequences on human 
health (Ware et al., 1986; US EPA, 2011) and on the ecosys- 
tem acidification (Schwartz, 1989). With a mean lifetime of 
few days in the troposphere (Lee et al., 201 1 ; He et al., 2012), 
emitted SO 2 is quickly oxidized to form sulfate aerosols. 
The resulting aerosols exert influences on the atmospheric 
radiative balance and cloud microphysics (e.g., McFiggans 
et al., 2006). SO 2 is emitted into the atmosphere mainly 
from anthropogenic sources such as fossil fuel combustion 
and industrial facilities. In the US these emissions repre- 
sent more than 90% of SO 2 released into the air (US EPA, 
2011). Since the implementation of national environmental 
regulations (e.g. 1990 Clean Air Act Amendments in the 
United States), a significant decrease of these emissions has 
been observed over the past 30 years. To keep track of SO 2 
emissions, this gas is monitored throughout the country by a 
system of continuously sampling ground-based instruments, 
and also by episodic intensive field campaigns. These cam- 
paigns are particularly valuable because the instruments de- 
ployed on the ground and from aircraft give not only the op- 
portunity to validate and improve the ability of space-based 
instruments to monitor air pollutants, but also provide the op- 
portunity to evaluate chemical transport models that simulate 
the SO 2 and sulfate lifecycle (Easter et al., 2004; Liu et al., 
2005; Goto et al., 201 1). The purpose of this paper is to take 
advantage of the data measured during the Frostburg field 
campaign held in Maryland during November 2010 to eval- 
uate the S0 2 simulated with the GEOS-5/GOCART model. 
We first describe in Section 2 the aerosol model and give a 
brief description of the SO 2 sources and the chemical pro- 
cesses considered within the model. In Section 3 we start 
by validating the modeled SO 2 at the surface over the con- 
tinental US using the data collected by EPA. In Section 4 
we evaluate the GEOS-5 simulated SO 2 with measurement 
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data taken during the campaign. Section 5 reports the con- 
clusions. 


2 Representation of Aerosols in the GEOS-5 Earth 
Modeling System 

The Goddard Earth Observing System version 5 (GEOS-5) 
model, the latest version from the NASA Global Modeling 
and Assimilation Office (GMAO), is a weather and climate 
capable model described by Rienecker et al. (2008). The 
GEOS-5 system includes atmospheric circulation and com- 
position, oceanic and land components. By including an 
aerosol transport module based on the Goddard Chemistry 
Aerosol Radiation and Transport (GOCART) model (Chin et 
al., 2002), GEOS-5 provides the capability of studying atmo- 
spheric composition and aerosol-chemistry-climate interac- 
tion (Colarco et al., 2010). In addition to providing reanaly- 
ses of traditional meteorological parameters (winds, pressure 
and temperature fields, Rienecker et al. (2008)), the inclusion 
of aerosols provides the background information for GEOS-5 
to produce reanalyses of aerosol fields using retrieved aerosol 
optical depth (AOD) from the space-based instrument Mod- 
erate Resolution Imaging Spectroradiometer (MODIS). The 
GEOS-5 near-real time system runs at a nominal 25 km hor- 
izontal resolution with 72 vertical levels between the surface 125 
and about 80 km. For this study, the model was run at various 
horizontal resolutions, 0.25° x 0.315° with sensitivity exper- 
iments also carried out at 0.5° x 0.625° latitude by longitude. 

GEOS-5 can be run in climate simulation, data assimila- 
tion, or replay modes. In the data assimilation mode, a mete- 130 
orological analysis is performed every six hours to constrain 
the meteorological state of the model. In the replay mode, a 
previous analysis, generated with the same version of model, 
is used to adjust the model’s meteorological state much like 
a Chemical Transport Model (CTM) with the difference that 135 
in GEOS-5 the aerosol transport dynamics are entirely con- 
sistent with the model thermodynamical state at every time 
step between analysis updates. For our replay simulations 
the GMAO atmospheric analyses from the Modern Era Ret- 
rospective analysis for Research and Applications (MERRA) 140 
(Rienecker et al., 2011) available every six hours are used. 

The GOCART module simulates five aerosol types: dust, 
sea salt, black carbon, organic carbon and sulfate aerosol. 
The sulfur chemistry processes considered are based on Chin 
et al. (2000a). Sulfate aerosol is mostly formed from the 145 
oxidation of SO2. All simulations include emissions of 
dimethysulfide (DMS), SO2 and sulfate and we use pre- 
scribed oxidant fields (hydroxyl radical (OH), nitrate radi- 
cal (NO3) and hydrogen peroxide (H2O2)) from a monthly 
varying climatology produced from simulations in the NASA 150 
Global Modeling Initiative (GMI) model (Duncan et al., 
2007; Strahan and Douglas, 2004). A small amount of SO2 
is produced by the oxidation of DMS, which is emitted nat- 
urally from marine phytoplankton. We use a monthly vary- 
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Fig. 1. SO 2 emissions released by coal fired plants in 
2007 over the United States (EPA source available at 
http://www.epa.gov/cleanenergy/energy-resources/egrid/). The 
circle size is proportional to the emission rates. 


ing climatology of oceanic DMS concentrations (Kettle et 
al., 1999), with emissions calculated using the surface wind- 
speed dependent (Liss and Merlivat, 1986) parameterizations 
of air-ocean exchange processes. The main source of SO2 
is anthropogenic, mainly from fossil fuel combustion from 
power plants and industrial activities (US EPA, 2011). 

Figure 1 maps the emissions of SO2 released from coal 
fired power plants (in tons) over the US in 2007. In this study, 
two different data sets of anthropogenic emissions and two 
assumptions about the injection height are considered in our 
simulations to assess the effect of the emissions on SO2 sur- 
face concentration. At the time of the campaign, the annual 
anthropogenic emissions of SO2 were taken from Streets et 
al. (2009). In the GEOS-5 control simulation, this emission 
was injected into the lowest model level. All simulated re- 
sults using this configuration are hereafter called the ’’Con- 
trol Run” or CR. 

Recently, a new Emission Database for Global Atmo- 
spheric Research (EDGAR) version v4.1 dataset (European 
Commission, 2010) became available at 0.5° horizontal res- 
olution and has the advantage of providing the 2005 anthro- 
pogenic emissions of SO2 by source categories. This new 
set of emissions allowed us to emit the non-energy emissions 
(from transportation, manufacturing industries, residential) 
into the lowest GEOS-5 layer and the energy emissions from 
power plants at higher levels between 100 and 500 meters 
(between the 2 nd and 4 th model layers). The results are 
herein referred to as the ’’Revised Run” or RR. 

Figure 2 shows a comparison of the SO2 anthropogenic 
emissions by source category: energy- source sector and non- 
energy- source sector, based on the EDGAR 2005 database as 
used in our revised simulation. Most SO2 emissions are re- 
leased from power plants, so it is important to consider the 
emission injection above 100 m due to the stack height and 
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Fig. 2. SO 2 anthropogenic emissions from the EDGAR v4.1 regrid- 
ded at 0.25°x0.3125° resolution in 2005 for non energy and energy 
sectors. 
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Fig. 3. GEOS-5/GOCART monthly mean of SO 2 surface-level (re- 
vised run) for January and July 2010. 


plume rise. We assume these emissions are constant through- 
out the year. Furthermore, other anthropogenic emissions 
include aircraft and ship traffic emissions from Mortlock et 
al. (1998) and Eyring et al. (2005) respectively. We assume 
3% of the SO 2 anthropogenic emissions are directly emitted 
as sulfate. All the simulations include also biomass burning 
emissions of SO 2 following the Quick Fire Emission Dataset 
(QFED) inventory and SO 2 emissions from continuously 
eruptive volcanoes that are based on data from the Global 
Volcanism Program database (Siebert et al., 2002) and To- 
tal Ozone Mapping Spectrometer (TOMS) and Ozone Mon- 
itoring Instrument (OMI)’s SO 2 retrievals (Carn et al., 2003; 
Krotkov et al., 2006) while emissions from explosive volca- 
noes follow the Aerocom inventories (Dentener et al., 2006). 
SO 2 is removed in the atmosphere by dry and wet deposition 
and oxidized to sulfate by chemical reaction. The main ox- 
idation pathways are in the aqueous phase by H 2 O 2 and in 
the gas phase by OH (Chin et al., 2000a). We save the model 
tracer fields every three hours during our simulation. Figure 3 iss 
shows results of the simulated SO 2 surface concentrations for 
January and July 2010. The highest SO 2 concentrations are 
found over eastern Asia, Europe, and North America, which 
are major anthropogenic source regions. SO 2 concentrations 
are higher during the winter; this seasonal variation can be 190 
explained by the seasonal SO 2 oxidation rates, which are 
slower in winter than in the summer (Chin et al., 2000b). The 
planetary boundary layer (PBL) dynamics is also responsible 
for this seasonal cycle of S02 concentrations. Figure 4 shows 
an evaluation of the GEOS-5 simulation of the SO 2 lifetime 195 
in black by comparison with the analysis made by Lee et 



Fig. 4. GEOS-5/GOCART monthly SO 2 lifetime for the year 2010 
compared to the study made by Lee et al. (2011) over the eastern 
United States (35.2°N - 44.5°N, 68.4°W - 81.6°W) during daytime. 


al. (2011) with the GEOS-Chem chemical transport model 
in red and in- situ measurements-based lifetime in blue. The 
mean S02 lifetime from GEOS-5 simulations are calculated 
over the eastern US (35.2°N - 44.5°N, 68.4°W - 81.6°W) 
and during daytime as Lee et al. (2011) but for the year 2010. 
The seasonal variation of the SO 2 lifetime from GEOS-5 is 
globally consistent with the seasonal variation found with the 
GEOS-chem model and the in-situ measurements. While the 
mean SO 2 lifetime from GEOS-chem are generally shorter 
than the in-situ measurement-based lifetime, the mean SO 2 
lifetime from GEOS-5 simulations are generally higher than 
the in-situ measurements, except during the winter. However, 
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the GEOS-5 SO 2 lifetime values are quite close or within the 
range defined by the uncertainty interval of in-situ measure- 
ments. The differences in the transport and in the emissions 
are among the possible reasons that may explain the discrep- 
ancy with the GEOS-Chem model. In addition the oxidant 
fields in GEOS-5 are not interactive and depend instead on 
fields from a different model from a different period. 


3 Model comparison to EPA surface measurements 

In this section we evaluate the modeled surface concentra- 
tions of SO 2 and sulfate over the US for the control and re- 
vised runs for the year 2010. For this study we used data col- 
lected by EPA, local and state control agencies which main- 
tain air quality monitoring networks over the US available 
from the EPA Air Quality System (AQS) (US EPA, 2010). 

3.1 Sulfur dioxide 

Figure 5 shows the SO 2 daily mean comparisons for the con- 
trol run (top) and the revised run (middle). The ”EPA” daily 
averages of SO 2 concentration were calculated using hourly 
concentrations collected from 102 sites obtained from the 
EPA AQS. A kernel density estimation (KDE) (Silverman, 
1986; Scott, 1992) was applied to approximate the joint prob- 
ability density function (PDF) of observed and modeled SO 2 
daily mean surface concentrations. Since SO 2 is usually log- 
normally distributed, the correlation coefficient (r), the Root 
Mean Square of the differences (GEOS-5-EPA) (RMS), the 
standard deviation (STDV) and the mean differences are cal- 
culated for logarithmically transformed data (summarized 
in Table 1 as well as the parameters in the original units 
calculated using the equations described in Limpert et al. 
(2001) (Appendix A)). For both plots, the scatter between 
modeled and observed daily means is significant with cor- 
relation coefficients, r=0.49 and r=0.42 for the control and 
revised run respectively. However, the agreement between 
the observed and modeled daily mean is better with the re- 
vised run, with lower values for the RMS and the mean dif- 
ference. The STDV is almost the same for both the con- 
trol and revised runs. One of the reasons for this discrep- 
ancy might be attributed to the change in absolute magni- 
tude of the SO 2 emissions datasets used in the control and 
revised runs, but we noticed only small differences between 
the two datasets. Another plausible explanation is the emis- 
sion injection height considered in the model. The verti- 
cal placement of emissions in the revised run decreases the 
high bias between observations and simulations at the sur- 
face. The remaining bias between observations and revised 
model SO 2 simulations may be explained by the error of rep- 
resentativeness associated with the incompatibility between 
in-situ measurements and grid-box mean values predicted by 
the model. As an attempt to filter out the in-situ measure- 
ments that are very unrepresentative of the grid-box mean 250 


conditions, the bottom plot of Fig. 5 presents the results af- 
ter a statistical quality control was performed with the adap- 
tive buddy check of Dee et al. (2001). For a given ob- 


Fig. 5. Comparison of daily averaged surface SO 2 concentration in 
2010 for 102 EPA sites. The model results are from the control run 
(top), the revised run (middle) and after the adaptive buddy check of 
Dee et al. (2001) was performed on the model revised simulations 
(bottom) (RR/bc). 

servation, this method consists of looking at nearby model- 
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Table 1. Summary of SO 2 surface comparison results (n is the number of points; r is the correlation coefficient; STDV is the standard 
deviation and diff is the mean difference in the logarithmic scale, the parameters with a ’*’ are the values in the original data scale as 
described in Limpert et al. (2001), Appendix B). 



n 

r 

STDV 

STDV* 

RMS 

RMS* 

diff 

diff' 



(log) 

(log) 

(ppb) 

(log) 

(ppb) 

(log) 

(ppb) 

Control run 

24916 

0.49 

1.00 

5.64 

1.38 

7.08 

0.95 

4.29 

Revised run 

27435 

0.42 

1.04 

3.22 

1.08 

3.95 

0.29 

2.29 

Revised run/buddy-check 

22538 

0.66 

0.73 

1.27 

0.74 

1.98 

0.15 

1.52 
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Fig. 6. The first column is r, the STDV and the absolute value of the mean difference between the modeled (control run) and observed daily 
averaged SO 2 surface concentrations for each SO 2 EPA site in 2010. The second column is the change in r, STDV, and absolute value of the 
mean difference for the revised run relative to the control run. The third column is the same, but showing the difference between the revised 
run (with buddy check of Dee et al. (2001)) and the control run. The color coding in the second and third column is such that blue indicates 
improvement relative to the control run. 


observations discrepancies and discarding those observations 
that cannot be corroborated by their neighbors. A brief sum- 
mary of the algorithm is given in Appendix B. After remov- 
ing observations that failed this adaptive buddy check (Fig. 5 
255 - bottom plot), the new comparison is quite improved with 270 

r that increased and is equal to 0.66 and lower values of 
the RMS, SDTV and the mean difference. The explana- 
tion for the remaining bias observed after the quality con- 
trol could be the year (2005) of the emission dataset with 
260 emissions too high for the year 2010. According to EPA 275 
(e.g., http://www.epa.gov/air/airtrends/sulfur.html) the aver- 
age SO 2 concentrations have decreased substantially over 
the years because of the application of SO 2 control mea- 
sures. Based on 341 US monitor sites, a 60% decrease in 
265 national average was found between 2000 and 2010. If we 


look site by site, Fig. 6 presents the change in the r (top), the 
STDV (middle) and the absolute value of the mean difference 
(bottom) between modeled and observed daily averaged sur- 
face SO 2 for the control run on the left, the revised run in the 
middle and after the buddy check on the right. While the cor- 
relation coefficient increased from values lower than 0.4-0.6 
for the control run to values greater than 0.6 after the buddy 
check, we see that the STDV increased over New England 
and slightly decreased elsewhere for the revised run, the de- 
crease is more significant after the buddy check. Concerning 
the absolute value of mean difference, we notice a decrease 
more and more significant between the control, the revised 
run and after the buddy-check. 
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3.2 Sulfate aerosol 

Figure 7 shows comparisons similar to the ones on Fig. 5, 
but for sulfate. The daily means are directly provided by 
the EPA AQS and are available every one, three or six days 
for a total of 250 sites. Figure 7 includes also a compari- 
son with the sulfate simulated with the GEOS-5 aerosol as- 
similation system, assimilation of MODIS AOD in the re- 
vised version of the model has been performed. On average 
the modeled sulfate concentrations are higher than the ob- 
servations, regardless of the model or data assimilation sys- 
tem used. The values of r, the RMS, STDV and the average 
differences are slightly different for the control, revised sim- 
ulations and the reanalysis (summarized in Table 2). This 
suggests that the SO 2 emissions injections as well as the as- 
similation of AOD observations into the model have a low 
impact on the daily mean sulfate comparisons. Eike for the 
SO 2 study, the measurements have been quality controlled 
using the buddy-check scheme (Fig. 7), permitting an in- 
crease r from 0.71 to 0.79, the RMS, the STDV and the mean 
difference have been divided by almost a factor 2. Coupled 
with the longer lifetime of SO 2 in Fig. 4 and 5 and, hence, 
too slow production of sulfate, our results suggest we may 
strongly underestimating the losses of sulfate aerosol. When 
looking site by site (Fig. 8), while the values of r decrease 
with the revised simulations for some sites, the application of 
the buddy check lead generally to greater and significant cor- 
relation coefficient values; the STDV values have not really 
changed between the control and revised runs but the values 330 
tend to decrease after the buddy check. Finally we see also an 
improvement in the absolute values of the mean differences 
after the revised and more importantly after the buddy check 
simulations. 


4 Evaluation of SO 2 in the model: comparison with 
measurement data during the Frostburg campaign in 
Maryland 340 

In Section 4 we concentrate our evaluation of the model 
performance in a smaller region using data collected dur- 
ing an air quality campaign in western Maryland in Novem- 
ber 2010. The Frostburg campaign was a regional air qual- 345 
ity campaign conducted by investigators from Washington 
State University (WSU), the University of Maryland (UMD) 
and the NASA Goddard Space Flight Center (GSFC) dur- 
ing two weeks in November 2010. The campaign took place 
in Western Maryland and provided direct measurements of 350 
SO 2 among other atmospheric constituents. The interest of 
this region is based on the abundance of SO 2 from the Ohio 
River Valley, surrounded by several power plants (Figure 9). 

In this section, we make use of several data sets available 
during this campaign to evaluate the anthropogenic SO2 con- 355 
centration simulated by GEOS-5. 



Fig. 7. Comparison of daily averaged sulfate surface concentrations 
for 250 EPA sites in 2010. The model results are from the control 
run, the revised run, the aerosol assimilation system and the revised 
simulations combined with the buddy check of Dee et al. (2001). 


4.1 Surface analysis: comparisons at Piney Run Station 

The observed and simulated monthly mean SO 2 at the sur- 
face at Piney Run station are shown in Figure 10. This site 
is located in a mountain valley close to Frostburg, and is an 
ideal location for SO 2 monitoring due to its close proximity 
to power plants stations, with the nearest one, Warrior Run, 
located south of Cumberland. Globally, the model captures 
the observed month-to-month variability of SO 2 with a win- 
ter maximum for both the control run in red and the revised 
run in black, as stated in section 2, the oxidation rates and the 
PBF dynamics are responsible for this seasonal variation. 

In the control run (the red line in Figure 10), we see that 
the model overestimates the observed SO 2 values by a factor 
of 4-5. This result is consistent with the general findings of 
section 3: the revised vertical placement of SO 2 emissions 
has a positive impact on the simulated surface values of SO 2 . 
This is shown with the revised run (in black) where the model 
values are in better agreement with the observations and the 
overestimation is less than a factor 2. Fike seen previously, 
an explanation of the positive bias remaining might be at- 
tributed to the 2005 emissions inventory and the recent de- 
creasing trend of SO 2 pollution over the US noted by EPA. In 
particular in Piney Run, the concentrations of SO 2 decreased 
50% between 2006 and 2010. 

Figure 11 shows the comparison of the daily mean SO 2 
surface concentrations to the measurements at Piney Run 
during 2010. Again, we see the better agreement between 
the revised run and the observations. 
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Table 2. Summary of sulfate comparison results. 



n 

r 

STDV 

STDV* 

RMS 

RMS* 

diff 

diff' 



(log) 

(log) 

(Mg/m 3 ) 

(log) 

(Mg/m 3 ) 

(log) 

(Mg/m 3 ) 

Control run 

17707 

0.71 

0.70 

1.54 

0.81 

2.46 

0.41 

1.92 

Revised run 

19658 

0.64 

0.80 

1.95 

0.90 

2.84 

0.40 

2.06 

Revised run / aerosol reanalyses 

19657 

0.65 

0.79 

1.99 

0.91 

2.92 

0.44 

2.12 

Revised run / buddy-check 

16444 

0.79 

0.54 

0.81 

0.57 

1.62 

0.19 

1.40 
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Fig. 8. Same as figure 5 but for daily averaged sulfate surface concentrations. 


5 Column amount analysis: comparisons to a MF-375 

DOAS instrument 

Simulated SO 2 column amount is evaluated with measure- 
360 ments from the Multifunction Differential Optical Absorp- 380 
tion Spectroscopy (MFDOAS) instrument developed at WSU 
(Herman et al., 2009; Spinei et al., 2010), deployed on the 
roof of a building at Frostburg State University (FSU) for the 
campaign. This instrument measures the direct sun irradi- 
365 ance and scattered sunlight in spectral UV and visible wave- 385 
lengths 281 - 498 nm at 0.83 nm spectral resolution recorded 
simultaneously with a CCD detector in the spectrograph fo- 
cal plane. Analysis of the measured spectra is done using the 
DOAS technique which is based on the Beer-Fambert law 
370 which states that the relationship at a wavelength between 390 
the intensity of the incident solar light and the transmitted 
one attenuated due to absorption and scattering by aerosols 
and molecules in the atmosphere (e.g., Platt, 1994; Plane and 
Smith, 1995). SO 2 column density is measured with an un- 


certainty less than 0.03 DU. A description of this instrument 
as well as the DOAS technique can be found in Spinei et al. 
(2010). Figure 12 shows the comparison between the col- 
umn density measured by the MFDOAS and simulated by 
GEOS-5 during daylight hours from 13:30 UTC until 21:00 
UTC on November 08 and 09. We notice that changing from 
one emission dataset to the other shows not much change 
on the total column amount between the two runs; it con- 
firms the small changes in the absolute magnitude of the SO 2 
emissions between the two datasets. Accounting for the un- 
certainty on the ground-based instrument, the comparison is 
rather satisfying with both the control and revised run but 
we notice that the model does not reproduce the observed 
diurnal variations. Besides the lack of diurnal variation in 
the prescribed emissions, an explanation might be the spatial 
resolution of the model (^25 km) and the offset pointing of 
the MFDOAS instrument when looking at the sun. 
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Fig. 9. Frosburg campaign regional map.Yellow circles are coal 
fired power plant stations; the circle size is proportional to the emis- 
sion rates. Piney Run station denoted by the white symbol o is lo- 
cated at 39.70°N and -79.01°W. Cumberland (in pink) is located at 
39.62°N and -78.77°W and the Frostburg State University (in blue) 
is located at 39.65°N and -78.93° W. Flight track on 11/08/2010 is 
in red. 



20 



Fig. 11. Times series of daily averaged concentrations of SO 2 at the 
surface in 2010 at Piney Run station. Blue squares are observations, 
red circles are model simulations with the control run, black circles 
are revised model simulations. 


11/08/2010 



Fig. 12. Daily variations of the SO 2 total column amount on 
1 1/08/2010 (top) and 1 1/09/2010 (bottom) at the Piney Run Station. 
Blue squares are MFDOAS measurements, red and black circles are 
model simulations with the control and revised runs respectively. 


Fig. 10. Monthly averaged concentrations of SO 2 at the surface in 
2010 at the Piney Run station. Blue squares are observations, red 
circles are model simulations with the control run, black circles are 
revised model simulations. Vertical bars are the standard deviations 
of monthly values for the model, shaded blue area for observations. 

5.1 Vertical analysis: comparisons to aircraft measure- 
ments 

The GEOS-5 simulated vertical distribution of SO2 is com- 
395 pared to aircraft measurements conducted on two different 
days during the campaign. The flights were made on the 
UMD Cessna 402B aircraft, which was equipped with a mod- 415 
ified pulse-fluorescence instrument to measure the in situ 
SO2 concentration (Taubman et al., 2006). The aircraft flight 
400 path on November 8 is shown on Fig. 9. Important regional 
power plants are marked by yellow circles in Figure 9, with 


the size of the circle indicating the magnitude of SO2 emis- 
sions. November 8, 2010 featured sustained winds as high 
as 29 km/h with gusts to 45 km/h around the time of the 
405 flight. November 9, 2010 was considerably calmer, with 
sustained winds under 19 km/h and gusts noted over Cum- 
berland around the time of the flight. These information 
were recorded at the airport, which is not an official National 
Weather Service reporting station, but they were also backed 
410 up by the informal observations of the airplane’s crew. Both 
flights lasted about two hours and were characterized by 
spiraling climbs and descents over Frostburg (39.65°N, - 
78.93°W) and Cumberland, Maryland (39.62°N, -78.77°W). 
Figure 13 shows the simulated vertical profile of SO2 for the 
control (left) and revised (middle) runs sampled along the air- 
craft flight path, as well as the comparisons of the modeled 
SO2 concentration from the revised run only to the aircraft 
observations for both days. The dark black lines in Figure 
13 show the modeled SO2 extracted exactly at the aircraft 
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position, while the blue shading shows the range of the mod- 
eled SO2 concentrations for the surrounded grid boxes (25 
km in the horizontal direction and 200 meters in the vertical 
direction). If we look at the vertical profiles comparisons be- 
tween the control and revised runs, we notice small changes 
between the two runs. On November 8th, GEOS-5 captures 475 
most of the major features of the aircraft observations, in- 
cluding the sharp vertical gradient encountered as the aircraft 
made its vertical profile near Cumberland (at about 60 min- 
utes of flight time). The turbulent mixing and strong winds 
during this day explain the air well mixed and coming from a 
much larger area. On November 9th the model also captures 480 
many of the aircraft variations but misses the observed high 
values between 60-80 minutes flight time. During this time 
frame, the aircraft was flying over Cumberland, near the coal 
fired power plant Warrior Run. The calmer weather condi- 
tions during this day may explain the high values observed 
locally that could not be reproduced by the model with a 25 
km resolution. Concerning the simulated surface-level SO2, 
like seen in more details in sections 3.1 and 4.1 we notice a 
slight overestimation of the SO2 surface-level concentration 
at the beginning and at the end of the flight on both days. 435 

6 Conclusions 

The Frostburg campaign that took place in Maryland in 
November 2010 was a good opportunity to evaluate the SO2 490 
simulated by the GEOS-5/GOCART system. By comparing 
the modeled SO2 against observed data, such as aircraft and 
ground-based measurements from a ground-based system in 
Frostburg, we have first diagnosed that the SO2 concentra- 
tions was overestimated at the surface and adjusting the ver- 495 
tical placement of the SO2 anthropogenic emissions inside 
GEOS-5 improved the SO2 surface concentrations without 
changing considerably the integrated total column amount. 
The improvement in our treatment of the SO2 anthropogenic 
emissions was confirmed with the analysis performed over 500 
the US using the EPA ground-based measurements. 

The comparisons of the vertical profile with aircraft data 
showed that despite the spatial coarse resolution of GEOS-5, 
most of the major features of the aircraft observations were 
reproduced by the model on November 8 because the weather 505 
was dynamic with turbulent mixing and strong winds. In 
contrast the analysis on November 9 shows that during quiet 
days, GEOS-5 will have difficulty of detecting plumes, espe- 
cially in the vicinity of point source. Concerning the GEOS-5 
simulated sulfate, the comparisons with the EPA data show 5 io 
that the changes in the SO2 emissions dataset and vertical 
distribution did not affect much the simulation of the sulfate 
at the surface, the positive bias observed with the control run 
remains with the revised run. These comparisons suggests 
that there might have an underestimated loss of sulfate in the 
model. A full analysis of the chemical processes could not 515 
be performed with the available data and there is a possibil- 


ity that part of this process could also explain part of the bias 
remaining in the SO2 and sulfate comparisons. 


Appendix A 

The lognormal distribution 

A random variable X is lognormally distributed if Y = log(X) 
has a normal distribution. The mean X and the standard de- 
viation sx of the normal variable are related to the Y and sy 
of the lognormal variable by (Limpert et al., 2001) : 

X = exp(F + sy 2 /2) (Al) 


s x =X y / exp(s F 2 -l) (A2) 


Appendix B 
Adaptive Buddy Check 

In the buddy-check algorithm of Dee et al. (2001), first a 
background check is performed where differences between 
the observed and modeled daily means are analyzed in or- 
der to identify a set of suspect observations, given a specified 
tolerance. An iterative buddy-check is then performed on 
each suspect observation using the remaining reliable obser- 
vations (called ’’buddies”) within a specified radius to per- 
form a refined acceptance test. The tolerance used for this 
buddy check is adaptive in the sense that current values of 
the observation minus model departures are used as a local 
modulator of the innovation variances used in the thresh- 
old test. Notice that before applying the buddy check the 
observation-model departures must be unbiased by removing 
the mean value. Figure B 1 shows the PDF of the points re- 
moved after the buddy check is performed for SO2 . Although 
in some cases GEOS-5 simulates lower SO2 surface values 
than the ground-based measurements, the majority of points 
removed after the buddy check are due of an overestimation 
of the GEOS-5 simulations compared to EPA measurements. 
While misplacement of plumes by the model could account 
for some large discrepancies that would be flagged by the 
buddy check, there is no reason to expect that these discrep- 
ancies would be of a given sign. Therefore, the positive bias 
of the removed observations may point to excessive emis- 
sions by GEOS-5 at specific locations. 
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Fig. 13. GEOS-5/GOCART SO 2 simulations for the control run (left), the revised run (middle) along the flight track on 1 1/08/2010 (top), on 
11/09/2010 (bottom). Modeled SO 2 vertical profiles for the control run (left) and revised run (middle), the white line is the aircraft altitude, 
on the right, the red line is the observed SO 2 concentration, the black line is the modeled SO 2 concentration (revised run), and the blue 
shading shows the range of simulated SO 2 for the surrounded grids. 
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